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Abstract 

We use a mix of field theoretic and holographic techniques to elucidate various properties 
^ of quantum entanglement entropy. In (3 + l)-dimensional conformal field theory we study 

^ the divergent terms in the entropy when the entangling surface has a conical or a wedge 

^ singularity. In (2 + 1) -dimensional field theory with a mass gap we calculate, for an arbitrary 

^ smooth entanglement contour, the expansion of the entropy in inverse odd powers of the mass. 

^ We show that the shape-dependent coefficients that arise are even powers of the extrinsic 

^ curvature and its derivatives. A useful dual construction of a (2 + 1) -dimensional theory, 

which allows us to exhibit these properties, is provided by the CGLP background. This 
smooth warped throat solution of 11-dimensional supergravity describes renormalization 
group flow from a conformal field theory in the UV to a gapped one in the IR. For this flow 
we calculate the recently introduced renormalized entanglement entropy and confirm that it 
is a monotonic function. 



May 2012 

* On leave from Department of Physics and Center for Theoretical Science, Princeton University. 



Contents 



1 Introduction [T] 

2 The (2 + l)-dimensional entanglement entropy in free massive theories |4] 

3 Holographic computation of the (2 + l)-dimensional entanglement entropy 

in gapped backgrounds [7] 

3.1 IR behavior of the EE for a circle M 

3.2 IR behavior of the EE for a general entangling surface [ID] 

4 An example: CGLP background of M-theory [T2] 

4.1 The renormalized entanglement entropy in the UV and the IR [14] 

4.2 Tests of the shape dependence of the entanglement entropy [TH] 

5 Shape dependence of the entanglement entropy in (3 + l)-dimensional CFT [T71 

5.1 The entanglement entropy for a wedge [T7] 

5.1.1 The free scalar field [17] 

5.1.2 A holographic computation [19] 

5.2 The entanglement entropy for a cone [21] 

A Perturbed CFT ^ 



1 Introduction 

The ground state entanglement entropies have emerged as a useful set of quantities for 
probing quantum entanglement and the degrees of freedom of many-body ground states 
(see (l}]6| for reviews and references to earlier work). If we consider the entanglement entropy 
(EE) of a (i-dimensional spatial region and its complement, then the leading term is typically 
proportional to the area of the {d— l)-dimensional boundary in units of the lattice spacing e. 
The useful information is then encoded in the sub-leading terms which depend on the shape 
of the boundary. For example, in (3 + l)-dimensional CFT, it has been found [T] that the 
expansion of the entanglement entropy (EE) for a smooth closed entangling surface S has 
the simple geometrical structure. 
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where /c|Jj, = — 7^7^Vpn^ is the second fundamental form, 7^j, = g^^j, — n^^nl is the induced 
metric (first fundamental form) on S, and Rj] is the Ricci scalar of E, which equals twice its 
Gaussian curvature. The Weyl anomaly coefficients a and c are normalized above in such a 
way that (a, c) = {tiq + llni/2, no + 677-1/2), where no and ni/2 are the numbers of real scalar 
and Dirac fields, respectively. 

In (2 + 1) -dimensional CFT the structure of the entanglement entropy for a smooth 
contour E in a plane is 

S = J-j-F, (1.2) 

where £s is the length of the contour. There is no known expression for F in terms of the 
curvature of the boundary. However, if the boundary contains a cusp of length r^ax and 
opening angle f2, then S contains an additional singular term — /cusp(fi) log(rmax/e) [8-10 



In both field theoretic ^ and holographic [9, 10 calculations, /cusp(^) turns out to be a 
smooth convex function that interpolates monotonically between ~ l/i7 behavior at small 
angles and zero at Q = n. However, the details of the function are not universal — the 
holographic, free scalar, and free fermion calculations produce different functions /cusp(^)- 

In this paper we present new results on the shape dependence of entanglement entropy in 
two and three spatial dimensions, studying both smooth and singular boundary geometries. 
Many of our calculations rely on the geometrical approach to the calculation of entanglement 



entropy 11-14 based on the gauge/gravity duality 15 -17 , but we also present some purely 
field theoretic arguments. In three spatial dimensions we will consider EE for a conical 
entangling surface with opening angle Q and show that the calculation in AdS:^ produces a 
term ~ log^(rmax/e)- For a wedge of length L and opening angle fi, we will show that 
EE contains a divergent term ~ f^edge{^)L/e. Surprisingly, we find /wedgc(^) = /cusp(^) 
both in the free scalar field theory and in the holographic calculations. 

In two spatial dimensions, the entanglement entropy S{R) across a circle of radius R 
has been invoked recently in deriving general constraints on renormalization group flows in 
(2 + l)-dimensional field theory. If such a field theory is conformal, then 

S{R) = a—-F. (1.3) 



As established in 18,19 the subleading /^-independent term is related to the regulated 



Euclidean path integral Z of the CFT on the three-dimensional sphere S^: F = — log|Z|. 



For a field theory that flows from a UV CFT to an IR CFT, it was conjectured 19-22 



that FuY > -Fir. Recently, an ingenious proof of this F-theorem was constructed using the 
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disk entanglement entropy |^. An important ingredient in the proof is the "renormahzed 
entanglement entropy" introduced in 124 



F{R) = -S{R) + RS'{R) 



;i-4) 



which is a finite function for theories that are conformal in the UV. For a CFT this function 



takes the constant value F from (1.3). An important property of J^{R) is that in the limit 
of large R it approaches Fir 24 . Furthermore, J^'{R) = RS"{R). It was shown in |23] 



that for any Lorentz invariant field theory S"{R) < 0. This demonstrates that J^{R) is a 
non-increasing function and therefore proves the F-theorem. 



In 24 ,25 , some holographic calculations of the renormahzed entanglement entropy J^{R) 
were presented. Similarly, we will calculate J^{R) for the smooth Cvetic-Gibbons-Lu-Pope 
(CGLP) solution 26 of 11-dimensional supergravity, which is a warped product of M?'^ and 



8-dimensional Stenzel space 



27 



e^. This smooth warped throat is similar to the 



KS background of type IIB string theory |28j. The CGLP background describes RG fiow 
from the CFT3 dual to AdS4^ x V5 2 in the UV (for its different field theoretic descriptions, 
see 29, 30)), to a gapped theory in the IR. The masses of some of the bound states in the 



CGLP theory were calculated in [31] . Using the holographic approach to the entanglement 
entropy, we will confirm that J^{R) for the CGLP background is a monotonic function that 
approaches zero as ~ 1/R for large R. This function exhibits an interesting second-order 
phase transition at a special value of i?, where the bulk surface reaches the bottom of the 
throat and its topology changes. Transitions of this type have been observed in earlier 



holographic calculations 24 32] (see also 13 14,33 35]). 



Generally, for theories with a mass gap of order m, the large R expansion of the disk 
entanglement entropy is expected to have the form |4,36 38 



S{R) 



a h /3m{27rR) - 7 + 27r 



00 

E 

n=0 



C-l~2n 



mR) 



2n+l ' 



1.5) 



where 7 is the topological entangelement entropy 39 ,40] (in the simple cases we will consider, 

-2n-i g^j^g related to the 



7 vanishes). Following |4,36,38 we will show that the terms ~ {mRy 
anomaly terms in 2n + 3 spatial dimensions. 

In order to gain better insight into the structure of entanglement entropy for gapped 
theories, we will generalize from a circle of radius R to an arbitrary smooth contour S. In 
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this case, the general structure of the EE in gapped (2 + l)-dimensional theories is 



oo ~J] 

5. = af + /3m£,-7 + E^' (1.6) 

n=0 

where the coefficients c5]^_2„ are integrals of functions of the extrinsic curvature and its 



derivatives 37 . The expansion of S^, has only odd power of 1/m because on dimensional 
grounds these terms are multiplied by even powers of the extrinsic curvature and its deriva- 
tives. Since in any pure state, and in particular in the vacuum, the EE of a region is equal 



to that of its complement, we have the symmetry k ^ —n 37 . Generalizing the arguments 



of g 36 38 , we give a prescription for calculating the c_j^_2„ for massive free scalar and Dirac 
fields. Using a holographic description of large theories with a mass gap, we calculate the 
coefficients c^^ and c^g explicitly. We check the infrared expansion for the specific case of 
CGLP background. 

2 The (2 + l)-dimensional entanglement entropy in free 
massive theories 

In this section we show how to calculate the 1/m expansion of the entanglement entropy 
for massive free scalar and fermion fields in (2 + l)-dimensions. We will take the entangling 
surface to be a smooth, closed curve Si of length and extrinsic curvature k in the t = 
slice of flat M^.i. 

More generally, one could consider the case where the (2 + 1) -dimensional spacetime 
is described by a general manifold Ai. Using the replica trick one is then able to show 



that the entanglement entropy has the large mass expansion of the form given in (1.6) with 



/3 = —{uq + ni/2)/12 (see, for example, [5,41,42 ). The integers Uq and ni/2 denote the 
numbers of real scalar and Dirac fields, respectively, in (2 -|- l)-dimensions. The coefficients 
c^i_2n known explicitly in the case where Si has vanishing extrinsic curvature \E 



41 



42 



We will henceforth take Ai = M^'^ and allow the surface Si, which is taken to lie in the 
t = plane, to have a non-trivial extrinsic curvature. We want to determine the coefficients 
c^i_2n terms of integrals of functions of the extrinsic curvature and its derivatives. Our 
approach to the computation follows that of Casini and Huerta [4,36 ,38 , who showed how 



to compute the coefficients = C-i/R in the special case where Si is a circle of radius R. 

The calculation proceeds by considering a higher, even dimensional QFT consisting of 
free fields in M?'^ x T^, with k > 1 odd and the symmetric /c-torus of a large volume 



4 



Vo1(T'') = . In the following argument one can replace by an arbitrary scalable k- 
dimensional smooth manifold. We give the free fields a small mass M, which will act as 
an infrared regulator for the conformal anomaly. We want to calculate the entanglement 
entropy in this theory across the (1 + fc)-dimensional surface Si+fc = Si x T'^, which fills the 
fc-torus and is described by the smooth curve Si in the t = plane of M?'^. We may Fourier 
decompose the field modes in the compact directions to obtain an infinite tower of massive 
(2 + l)-dimensional fields, with masses 



27r^ ^ ^ 



M'+r-^] n^eZ. (2.1) 



i=l 



The entanglement entropy in the (2 + + l)-dimensional theory then becomes equal to the 
sum over (2 + l)-dimensional entanglement entropies for massive fields across the curve Si. 
Taking the large L limit, the spectrum of masses becomes continuous and we find 

^ST'W = /J-^i-fl 1) /'^p/-'4r"(yM^) ■ (2.2) 



where e is the UV cut-off. We now substitute the expansion of S^^^\m) given in (1.6) into 
(2.2). We see that the term in the expansion of S^^^\'m) which goes as 1/m'^ determines 



the logarithmic conformal anomaly term in S^^^^^^^\m). Turning this argument around, 
suppose the entropy of the (2n + 4)-dimensional theory has the anomaly term 



^C.W, =C„:rMMe), (2.3) 

log 



then we can immediately read off the coefficient c^\_ 



■2n- 



~Ea _ 7r(27r)"(2n-l)!! (2„+4) 



^-l-2n — Vol(T2"+l) "^^2^+2 • (2-^) 

The above formula is slightly modified for fermions. Dirac fermions in {2n + 4) dimensions 
are in a 2"'~^^-dimensional representation, which after dimensional reduction reduces to 2""*"^ 



(2 + l)-dimensional Dirac fermions. Thus, the right hand side of (2.4) should be divided by 
2"+i ^Qj. j)ii-ac fermions. As a corollary to this argument, we see that the absence of the loge 
terms in odd dimensional CFTs implies that the IR expansion of S^^^\m) contains only 
odd powers of 1/m, in agreement with the arguments in 37 . 



Let's see how this works explicitly when n = 0. The expression for s^^^^ is given in 
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(1.1). The Euler number x(^2) vanishes for S2 = Si x S^. The two normal vectors to S2 
are within M^'^, which we write with coordinates 



(2.5) 



where 6 has period 2tt. One of the normal vectors to S2 is timelike, = (1,0,0,0), where 
the fourth component is in the direction of the of length L, and its second fundamental 
form vanishes. Suppose Si is defined by a curve r = R{6). Then, the other normal vector 
is spacelike, = (0, r, —rR'{6),0)/ y^r^ + R''^{0), and this gives a second fundamental form 
with non-vanishing component 



26 



R\e) + 2R'\e) - R{e)R"{e) 

(i?2(^) + (i?'(0))2)3/2 



(2.6) 



where k{6) is the extrinsic curvature of the surface Si in the plane. It follows that both 



Ki^K" K^Ku become k'^{9). This leads to 



480 



{no + 3^1/2 



(is K 



(2.7) 



and 



5. 



(2+1), 



a 



"^(^0 + ^1/2)^E1 ^0 + 3^1/2 
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480 m 



ds K^ + 0{l/m^) 



(2. 



where we stress that ^1/2 the number of (2 + l)-dimensional Dirac fermions. 

In principle the calculation of the higher order corrections to the entanglement entropy 



in powers of 1/m would proceed analogously. For example, to calculate the coefficient c„3, 
which gives the order 1 / correction to the entropy, we would need to first calculate 
with S4 = Si X T^. We then expect 



(5+1) 



,(5+1) 
'S4 



Vo1(T3 



{Aq no + Ai/2 nf^^) ^ ds + {Bq Uq + B1/2 nf^^) j 



(2.9) 



for some coefficients (Aq, ^1/2) and {Bo, -B1/2), which should be functions of the 6-dimensional 
anomaly coefficients. We use the notation n^y2 stress that this counts the number of 
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(5 + 1) -dimensional Dirac fermions. This then leads to 



A ^^1/2 \ / , 4^Ar ^^1/2 \ / , (dnV 
no + ni/2 1 (b dsn + I no + ni/2 1 f "-^ ) 



. (2.10) 



This formula is consistent with the general arguments in 37 



3 Holographic computation of the (2 + l)-dimensional 
entanglement entropy in gapped backgrounds 

The (renormalized) entanglement entropy may be calculated holographically by following the 



usual procedure for holographic entanglement entropy 11-14 . Consider a ((i+l)-dimensional 



large N field theory with a D-dimensional gravitational dual. While we will ultimately be 
interested in (2 + l)-dimensional QFT, for now we keep the dimension d general. As in 14 , 
let the gravitational background have the Einstein-frame metric 

ds^jj = a{u)[du^ + (3{u)dx'^dx^] + gijdy^dy^ , dx^dx^ = —dt^ + dr"^ + r'^dVl?^_^ , (3.1) 

with a{u) > and P{u) > and i, j = d + 3, ■ ■ ■ , D. The compact, internal [D — d — 2)- 
dimensional manifold is taken to have a volume 

D 



V{u)= / n dy'^/d^g, (3.2) 



i=d+3 

which is a function of the holographic radial coordinate u. We assume that u has a minimal 
value Uo where a p-sphere in the internal manifold shrinks to zero size, resulting in V{uo) = 0. 
At uq we assume that all supergravity fields are regular, which implies a{uo) and /3(mo) are 
finite. The coordinate u ranges from infinity in the far UV to Uq in the far IR. Such geometries 
typically describe confining gauge theories. 

We further assume that the gravitational theory approaches a conformal fixed point in 
the UV (u = oo), and we work in coordinates where 



lim a{u) = auY , lim V{u) = Vuv , Piu) = exp ( — ^ '^^ ) + . . . , (3.3) 

u— >oo 11— >oo \ -^UV / 

where auv and Vuv are constants, and Luv is the radius of AdSd+2- 

We want to calculate the entanglement entropy in the QFT across a codimension two 
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spacelike surface T^d-i- The entanglement entropy 11 14 is calculated holographically by 



finding the {D — 2)-dimensional surface Sd_2, which approaches T^d-i at the boundary of 
the bulk manifold, is extended in the rest of the spatial dimensions, and minimizes the area 
functional 



Sj 



1 



AG 



(D) 
N 



d^'''a\lG\ 



,{D-2) 



ind 



(3.4) 



where G[^^'^'' is the induced metric on Tj£)_2- The entanglement entropy is then given by 
the functional S^. evaluated at the extremum. 

A case of particular interest is when the region T^d-i is the {d — l)-sphere of radius R. 
Writing the radial coordinate r as a function of the holographic coordinate m, the induced 
metric on S£)_2 is 



dsl = a{u)[{l + (3{u){dury)du'^ + l3{uy{u)dnl^-^] + g^jdy'dy^ 



(3.5) 



which gives the following expression for the area functional in terms of the unknown function 
r(u): 



S{R) 



Vol(5' 



d-i\ r°° 



dnr''-\n)g(u)^l + I3(u)(d^r) 



g(u) = a'''\u)P^''-^'^'\u)V{u) 



In general we need to first solve the Euler-Lagrange equation. 



(3.6) 



{d-l)r''-\u)g{u)^l + P{u){duTY 



du 



r'^ ^{u)g{u)/3{u){dur) 



(3.7) 



for the function r{u), then evaluate the area functional in (3.6) on the solution with a UV 



cut-off u < u\jy, then use (1.4) to construct the finite renormalized entanglement entropy. 



For non-trivial backgrounds this must be done numerically. To solve the equation of motion 



(3.7), we also need to specify the boundary conditions. There are two types of solutions with 



different topologies. 

One of them, which we will call the cylinder-type solution, terminates ai u = uq where 
the volume of the internal space becomes zero: V{uq) = 0. One can find the form of the 
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solutions r{u) for u near Uq by expanding (3.7) around u = Uq: 



r{u) = ro + 



d 



ro > 0. 



(3.^ 



The other type of solution, which we call the disk-type solution, has a tip at u = Umin > 
Uq, where the radius of the sphere becomes zero: r{umm] 
behave like 



0. For u near Wmin, the solutions 



dg{u. 



2/3(Mmin)fi''(Mmin) + 9{Umm) /3' {Ur_ 



-(w-w^i„)^/2 + 0((w-w^i„)3/2). (3.9) 



3.1 IR behavior of the EE for a circle 

We may obtain the IR asymptotic behavior of the entanglement entropy for Ei = through 
an analytic procedure, and in doing so we show that the renormalized entanglement entropy 
approaches zero in the IR from above like 1/R, where R is the radius of the . Note that 
in this section we restrict to the physical dimension d = 2. In the following section we 
generalize the computation by allowing for a general entangling surface Si. 

For now, we take Si = 5*^ of radius R. We assume that at large R the solutions to the 
Euler-Lagrange equations will be of the form r{u) = R + 6{u)/R, with 6{u) independent of 
R. Expanding the Euler-Lagrange equation in powers of l/R, we find the equation 



d 

du 



which may be integrated to obtain 



5{u) 



[g{u)P{u)5\u)] = g{u) 



(3.10) 



du' 



1 



g{u')P{u') 



du"g{u"). 



(3.11) 



Mo 



Expanding the area functional in (3.6) and using the equation of motion in (3.10) gives 

2%[i.^ + i?(f^ci.,(.)-^) 

(3.12) 



SiR) 



AG 



N 



e 

1 r°° 



dug{u)P{u)[5'{u)Y + 0{R 



where we used the boundary conditions 5(uoo) = and S'{uq) = for the solution, which 
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make the surface term vanish. The UV cut-off e defined by 



— = auvexp ) . (3.13) 



-^uv / 



To compare with (1.6), we set the mass m to unity and use the dimensionless radius R 



for convenience. We then see that we can make the identifications 



-1 /"°° du 



r du ( r , , , ,x 2 

/ f \nf \ I du g(u ) (t dsn . 



(3.14) 



Notice that the coefficient /3 is finite and independent of the UV cut-off e. To calculate 
the coefficients 0^3 we must consider a more general entangling surface. This is because 
dn/ds = for the circle. In the following section we generalize the above calculation to allow 
for a general, smooth entangling surface, and in doing so we calculate c^g. 



3.2 IR behavior of the EE for a general entanghng surface 

We would like to repeat the calculation in the previous section allowing for a general spacelike 
entangling surface Si. While we believe that the computation can be carried out in full 
generality, it is enough to restrict ourselves to a closed curve Si that is a boundary of a star- 
shaped domain]^ Such a curve can be parameterized using polar coordinates by a function 
Ra{9). We write the entangling surface as Ra{9) = AR{9), with R{9) a smooth function 
and A > 1. The IR limit corresponds to A large enough such that the extrinsic curvature is 
small, K\{9) = A^^k{9) ^ 1, along the entire curve. 
The induced metric on the bulk surface T,d-2 is now 

dsl = a{u) [{1 + /3{u){durf)du^ + P{u){r\u, 9) + {derY)d9'^ 

. . (3.15) 
+ 2l3{u){d0r){dur)d9 du] + Qijdy'dy^ , 

where the radial coordinate r{u, 9) is taken to be a function of the holographic coordinate u 
and the angular coordinate 9. We require that limtt_5.oo r{u, 9) = Ra{9). The area functional 

star-shaped domain is a set 5* C M" with the property that there exists a point xq ^ S such that the 
hne segments joining xq to all other points in S are contained in S. 
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for the entanglement entropy may be written as 
1 



2tt 



de / dug{u)^{r\u,e) + (9,r)2) {I + (5{u){durf) - I5{u){der)\dury 



Uo 



(3.16) 



with g{u) and f3{u) defined as before. 

We assume that in the IR the solutions to the Euler-Lagrange equation give 

riu, 9) = ARi9) + + 0(1/ A') , 



(3.17) 



with 6{u,6) order A''. We Substitute the ansatz in (3.17) into the area functional in (3.16) 



and expand in inverse powers of A up to and including terms of order 1/A: 
1 



4G^^^ Jo 



2lT I'OO 

dO I du 

Juo 



Ag{u)^yR{ey + R'{ey 



g{u)R{e) f R{e)l5{u) 
A \^/RH0)Tm9) 



-.{djf + 2K{d)5{u,d) 



(3.18) 



+ 0(1/A3 



where the extrinsic curvature of the entangling surface, fi;(6'), is given explicitly in (2.6) 
Applying the variational principle to find the Euler-Lagrange equation for 5(u, 6) gives 



— [g[u)(5[u)du5{u,e)\ = ^ ^^"^ 



(3.19) 



which may be integrated to give 



K{e)y/R{e)^ + R'{ey 



2 roo 



R{e) 



1 



du - , , ^, , 



du"giu"). (3.20) 



We want to calculate the terms in the expansion of of order 1/A^. These terms are 



completely determined by the expansion of r{u,6) in (3.17) through order 1/A. Expanding 
the area function in (3.16) through order 1/A^ and evaluating on the solution for S{u,6) 



given in (3.20) allows us to determine the 0^3 coefficients in (1.6): 



c^Q = a 



(1) 



ds 



.(2) 



ds 



(3.21) 
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with 



,(2) 



du 



ds 



1 



H 2 



du g{u') 



32G 



■uo 

du 



du' 



du" g{u") , 



uo 



g{ufP{u) 



du' g{u'] 



UQ 



(3.22) 



4 An example: CGLP background of M-theory 



The CGLP background 26 of M-theory is the gravitational dual of a gapped (2 + 1)- 



dimensional field theory, which nicely illustrates the general features discussed in the previous 
section. The supergravity background is a warped product of M^'^ and an eight-dimensional 



Stenzel space 27 



E 

i=l 



2 2 

zf = e . 



(4.1) 



where e is a real deformation parameter. When e = this equation describes an eight- 
dimensional cone whose base is the Stiefel manifold V52. 



As explained in 26,43 , the Stenzel space (4.1) can be parameterized by a radial coor 



dinate r ranging from to 00 and the seven angles in ¥5^2- At r = a 3-sphere shrinks to 
zero size, and the r = section is a round 5"^. 



The 11-dimensional metric is of the form of the metric in (3.1) if we identify the holo- 
graphic coordinate u with r, where tq = 0, anc|^ 



air 



4 



c^{t)H{t) 



V{t) = -3'/'71^H'/%t){2 + coshr)3/« sinh^/^ (-) sinh^/^r) . 
2 \ 2 / 



(4.2) 



The functions H{t) and c(r) are defined by 



^ ' 8l7r4 



dt 



(2+coshT) 



1/4 {t^ - 1)5/2 



3^/^ cosh^ f 
^(2 + coshr)3/4 



(4.3) 



where is the number of units of asymptotic G4 flux. In particular, notice that V{t = 0) = 
0, which is a result of the vanishing 3-sphere. For more details on the CGLP background 



see, for example, 26 43 



We follow the conventions of 43 and work in units where e 
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We begin by studying the entanglement entropy and renormalized entanglement entropy 
in the simpler case where the entangling region Si is taken to be a circle of radius R. In this 



case the Euler-Lagrange equation for the function r(r) in (3.7) may be solved numerically 
with the boundary condition r(r = oo) = R. In practice, we cut the space off at some 
large tuv- For each R > there exists a value T^in{R), which is the smallest value of r 
for which the function r(r) is defined. There exists a critical value -Rent ~ -73 for which 
'"(''"min = 0) = 0. For R < Rcrit the solutions to the equation of motion describe surfaces of 



disk type that behave as in (3.9) for The topology of these surfaces is that of a 



disk times V^^2- The solutions for R > i?crit are surfaces of cylindrical type that stretch to 
the bottom of the Stenzel space and behave as in (3.8) for r near tq = 0. The topology of 



these surfaces is that of a circle times the Stenzel space. 

In Figure [l] (a) we plot the numerical solutions to the equation of motion for a range 



r(T) 
0.8 r 



(b) 



0.6 
0.4 
0.2 
0.0 
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(c) 
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Figure 1: (a) Numerical solutions to the equation of motion for the holographic entangling 
surface, given by r(r), in the CGLP theory. The dotted red line indicates the critical value 
i?crit, where the solutions change from disk-type to cylinder-type, (b) A zoomed-in plot of 



the UV region, with disk-type solutions, where we plot the AdS approximation in (4.4) in 



dotted black, (c) A zoomed-in plot of the IR region, with cylinder-type solutions, with the 



analytic approximation given by 6{t) in (3.11) plotted in dotted black 



of < i^crit and R > Rent- III the far UV the solution for r(r) should approach the AdS 
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solution 



r(r) = Vi?2 _ 25/23i/2e-3r/2 _ (4.4) 
In Figure [l] (b) we zoom in on some of the disk-type solutions in the far UV and plot the 



solutions along with the asymptotic in (4.4). In the far IR region the cylinder-type solutions 
should be well approximated by the function 5{t) in (3.11). In Figure [l](c) we plot some of 
the cylinder-type solutions along with the analytic approximation. 

As was discussed in section |3| to calculate the renormalized entanglement entropy it is 
sufficient to evaluate the entanglement entropy with a strict UV cut-off. We cut off the space 
at a large r value ruy. We then numerically integrate the area functional and differentiate 
it to construct J-". A plot of the renormalized entanglement entropy along the RG flow is 
given in Figure [2j 




1.5 



2.0 



2.5 



R 



3.0 



Figure 2: The renormalized entanglement entropy J^{R) along the RG flow in the CGLP 
theory plotted in orange. The left dotted black curve is the asymptotic UV approximation 



to J^{R) given in (4.7). The right dotted black curve is the IR approximation to J^{R) given 
in (4.8). The dotted red line marks the value R = Rent- 



4.1 The renormalized entanglement entropy in the UV and the 
IR 

In the far UV we can treat the CGLP M-theory background as a perturbation of the 



AdS^ X V5,2 background. From (A. 7), we know that the UV fixed point has a renormal- 
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ized entanglement entropy 



uv 



167rAr3/2 
27 



(4.5) 



where we have used Vol(V5,2) = 277r^/128. To describe the RG flow in the vicinity of the UV 



fixed point, it is convenient to use the effective 4-dimensional metric in the form of (A.l) 
A straightforward calculation shows that at small y the function f[y) has the expansion 



/(2/) = 1 + 2^332/3^4/3 + ...^ 



(4.6) 



which, using (A. 11), implies the RG flow is driven by an operator in the UV field theory of 



dimension A = 7/3, which is consistent with 43 . Using (A. 12), we then see that at small R 

167riV3/2 / 3 



HR) 



27 



1 - ^2^/332/3^4/3 + . . . 
7 



(4.7) 



This function is plotted together with the numerical solution in Figure [2j Note that it is 
a very good approximation to the actual renormalized entanglement entropy for R < Rent- 
We also see explicitly that djiJ-' = at R = 0. 



In the IR we may use (3.14) and (3.21) to get an asymptotic expression for the renor- 



malized entanglement entropy, which gives 



J" 



167rAr3/2 /0.1959 1.845 x lO'^ , 

+ + 0(1/ R^) 



27 



R 



i?3 



(4.8) 



This function is plotted in Figure [2| which shows that it is a good approximation to the 
actual renormalized entanglement entropy at large R. 



4.2 Tests of the shape dependence of the entanglement entropy 

In this section we will consider a more general spacelike entangling surface Si, specified by the 



function R{6) in polar coordinates. We want to check (3.20), which gives an approximation 



to the cylinder-type solutions in the far IR. In particular, this equation claims that the 
variation of the bulk entangling surface S2 away from the straight cylinder is proportional 
to the combination K{e){y/R{ey + R'{e)yR{e)). 

As an example, we consider the entangling surface Ei plotted in Figure [3} which has a 
small extrinsic curvature along the entire curve. A good way of measuring the accuracy of the 



analytic approximation in (3.20) is to compare the function Ruy{9) — Rm{6) as computed 



15 



160 




Figure 3: (a) A plot of an example entangling surface Ei, described by the function R\jy(9) 
in polar coordinates, in the CGLP theory, (b) The extrinsic curvature kuy{6) for the 
entangling surface Ruy{6). The extrinsic curvature is small over the whole curve. 

0.004 r 



0.002 



-0.002 
-0.004 
-0.006 



•.3 



V 



Figure 4: The function Ruy{9) — Rib.{9) in the CGLP theory, with R\jy(9) plotted in fig. [sj 
The solid orange curve is computed by numerically solving the equation of motion for the 
holographic entangling surface S2. The dotted black curve is an analytic approximation, 
which is equal to — 5(r = 0,6'), with 5(0,6') given in (3.20). 



both numerically and using (3.20), which is done in Figure |4j Here the function Ri^{6) 
is the profile of the cylinder when it reaches r = 0: -Rir(^) = r(r = 0,6'). The analytic 
approximation simply gives Ruy{9) — Rm{0) ^ — 5(r = 0,^^), with S{0,9) given in (3.20). 
The two curves match extremely well. 



16 



5 Shape dependence of the entanglement entropy in 
(3 + l)-dimensional CFT 



The entanglement entropy for a smooth entanghng surface S in a (3 + l)-dimensional CFT 



is given in (1.1). However, if the surface is not smooth, for example if it has conical or wedge 
singularities, then there may be additional contributions to ( |1.1 ). In this section we consider 
the entanglement entropy for a wedge and a cone in a (3 + l)-dimensional CFT through 
both field theoretical and holographic computations. We find that wedge entanglement 
entropy acquires a 1/e divergence not present in (1.1). The cone entanglement entropy has 
a log^ e divergence as predicted by (1.1), but its correct coefficient is twice smaller than for 



a regulated version of (1.1). 



5.1 The entanglement entropy for a wedge 

The wedge is the surface in (3 + l)-dimensions given by {x^,x^,x^) = (r sin0, r cos0, 2;), 
where < r < 00, = and Q, and z z + L. We have compactified the z direction on 
a large circle of length L to avoid unwanted infrared divergences. We begin by using the 
replica trick to calculate the entanglement entropy with this geometry for a free scalar field, 
and this is followed by a holographic computation. 



5.1.1 The free scalar field 

Using the replica trick one can show that the entanglement entropy for the massive scalar 



field is given by 44 



a-l 

S=- J2^ogZ, 

k=0 



(5.1) 



where is the partition function of a scalar field 0^ on 4-dimensional Euclidean space with 
boundary conditions 



X e A 



(5.2) 
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where A is the region bounded by the wedge. Since the theory we consider is free, the 
partition function Zk is obtained from the integral of the Green's function: 

^ logZfc = -l [ d^'+'XGkiXj) . (5.3) 



The Green's function is the two-point correlation function of the free massive scalar and is 
subject to the following conditions: 

(_A^ + m^) Gk{X, X') = 6{X - X') , 

^ ^ (5.4) 
lim Gfc((f,e),X') =e2"*^ lim Gfc((X,e),X') , xeA. 

e—^0+ €—^0- 

We expand the Green's function in Fourier modes along the z-direction. The problem of 
finding the 4-dimensional Green's function then reduces to that of finding the 3-dimensional 
Green's functions for a cusp entangling surface for a tower of massive fields, with masses 

M' = m'+(^]\ neZ. (5.5) 



Using the result in [8] for the Green's function with a cusp entangling surface in 3-dimensions, 
we find 

G,{X,X') = -J2J2 I dX^^^^,9n{x)9:{x) , (5.6) 



U 71= — OO 



where we use x to denote the 3-dimensional coordinates {t, x^,x'^). The (?„ are the eigenfunc- 
tions of the 3-dimensional Laplace operator (— + M^), whose eigenvalues we denote by 
(A^ + M^). Using spherical coordinates with {t,xi,X2) = {p cos 9, p sin 9 sin (p, p sin 9 cos (p), 
the eigenfunctions are given explicitly by 

Ji+J>^p) 

where J is the Bessel function of the first kind. The functions ipv are the eigenfunctions of 
the angular laplacian on the two-sphere. 



^nMOA) = + / d9d(t> sin{9)\ij,{9A)? = 1, (5- 
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subject to the boundary condition 



lim 'ip^{- + e,> 



e^'^'- lim 



0+ 



e [0, 



(5.9) 



Preforming the integral over A in (5.6) and taking a derivative of the partition function 
in (5.3) with respect to gives 45 



^ lo Z — ^ coth ( 
dm? ^ Am \ 2 



(5.10) 



The sum over the eigenvalues v is divergent and needs regularization. The computation of 
the sum was carried out in detail in [s] , and one finds that the regularized sum only depends 



on k and the angle of the cusp Vt. After integrating (5.10) with respect to m , we find that 



the entanglement entropy for a wedge has the angle dependent UV divergence 



wedge 



OL-1 



a — 



k=0 



(5.11) 



where the function fc^sp^'^\^) is the same function as for the cusp in (2 + l)-dimensional 
CFT It behaves as fc^sp^^\^) ~ 1/^^ when the angle is very small, while it becomes 
zero aX Q = TT where there is no cusp in the entangling surface. The function /cusp(^) is not 
universal and depends on the type of matter. 



5.1.2 A holographic computation 

Next we compute the holographic entanglement entropy for the wedge. To this end, we use 
the following AdS^ metric. 



ds^' 



dy"^ - dt^ + dr"^ + r'^dcf + dz^ 



r 



(5.12) 



For simplicity, we have set the AdS radius to 1. The central charges a and c of the dual CFT4, 



normalized as in (1.1), are then determined by the 5-dimensional Newton constant |46 

457r 



a = 3 c 



^(5) 



(5.13) 



z ~ z + L} at the AdS boundary 



The wedge is defined byS = {0<r< rmax, 4 
y = 0. The large radius cut-off rmax and the length L are introduced to regularize the volume 
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of the wedge. As usual, to introduce a UV cut-off we will restrict y > e. The entanglement 
entropy functional is given by 



S 



4GP 



dr 



^— 2 + r 2 {dry) ^ + {d^y) ^ 



(5.14) 



where we take the holographic coordinate ?/ to be a function of (r, 0). 

We must find the function ?/(r, 0) which minimizes the entanglement entropy functional 
and approaches the wedge at the boundary of AdS^. The scaling symmetry of the spacetime 
and the wedge implies the following ansatz for the minimal surface [9 10 : 



9{4>) 



With this ansatz the initial value problem becomes first ordeij^ 



(5.15) 



£/W(i + r 



go = giO), ^7'(0)=0, 



(5.16) 



It follows that the angle of the wedge determines 5^0 as 



2" 



dg- 



go 



(5.17) 



Integrating this equation we find that, as in the (2 + l)-dimensional cusp calculation [9, 10 
the limiting value where gQ = is Q = it . 



The entanglement entropy is then found by evaluating the regularized functional in (5.14) 
on the solution to the equation of motion: 



S 



2L 



r-max rr/e 

/ dgh{g,go) 



^'-j'TV soe go 



2L 



r 



- gl 



+ 



dg{h{g,go) ~ g) 



go 



{5.U 



AG 



(5) 
N 



^-/S>)7 + 0(.») 



^Note that one must first find tlie equation of motion for y(r, (j)) and tfien subsequently substitute tlie 



ansatz in (5.15) 
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The functions h{g,go) and /^tdge(^) defined by 



h{g,go) 



9\l + 9' 



f°° dr r 
9o- dg{h{g,go) - g) . 



90 



(5.19) 



90 



Just like the free scalar field calculation for the wedge, (5.11 ), the holographic result has 



an L/e divergence that was absent for smooth entangling surfaces]^ A numerical plot of 
the function /^tdge(^) shown in Figure [5| where it can be seen that it goes to zero as fl 



approaches tt. When Q is small /wtdKe(^) diverges as 



.(hoi) 
J wedge 



0.646 



(5.20) 



A pole at = also appeared in the field theory computation ( 5.11[ ). With that said, the 
function /^edge(^) different from that of the scalar field theory, /cusp'^'^''(^^), which appeared 



in (5.11) 



A surprising result, however, is that after an overall rescaling the function /^tdge(^) agrees 
with function /iu°p''(f2) in 10 describing the holographic cusp anomaly in (2 + l)-dimensions. 
The normalization of the function /wedge (^) depends on the choice of the UV cut-off scale e. 
We introduce the normalized function /wtdge(^) = /wedge (^) tuning the constant a such 



that /wtdge(^) agrees with /cusp (f^) in the limit of f2 — )■ 0. We find a ~ 1.11 numerically, and 



.(hoi) 



in Figure |5| we plot both /wedge(^) and /cusp (^^). The plot shows that in fact the normalized 



Khol) 



function /wtdgc(^) same (within the numerical accuracy) as /cusp (^^), although the 



(hoi). 



definitions (5.17) and (5.19) appear quite different from those for the cusp in 10 . For a 



free scalar field the function fl^cd^\^) ^"^^ wedge also turned out to be the same as 
ftp''\n)hr the cusp (see (|5.11[)). It is very interesting that the appropriately normalized 



functions f{^) agree for the cusp and wedge geometries both in the free and in the strongly 
coupled theories that we have studied. 



5.2 The entanglement entropy for a cone 

In this section we show that when the entangling surface in (3 + l)-dimensional CFT has 
a conical singularity, the entanglement entropy acquires a log^(rmax/e) divergence. We take 

^ Since on the sides of the wedge the extrinsic curvature vanishes, it is reasonable to think of this term 
in the entanglement entropy as due to the wedge singularity. 
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7r 



6: 
5 : 
4: 
3: 
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1 : 



0.5 1.0 1.5 2.0 



Figure 5: The entanglement entropy for the wedge has a 1/e divergent term, whose coefficient 
depends on the angle of the wedge. This coefficient /wtdge(^) given explicitly in (5.17) 



and (5.19), and its normalization depends on the UV cut-off e. The black dotted line is the 
normalized function /^tdgc(^) ~ '^/wedge(^) ^"^^ wedge with a ~ 1.11 and the orange line 
is the function /iusp (^^) for the cusp in (2 + l)-dimensions. 

the entangling surface to be the cone defined by (r, 6', 0) = (r, where < r < Tmax, 

~ + 27r is the azimuthal angle, and Q is the opening angle of the cone. The large distance 
cut-off Tmax rcgulatcs the area of the cone. 



To begin, we will evaluate (1.1) for this surface. Even though this equation is only 
valid for smooth entangling surfaces, it does provide a quick way of seeing how the log^ e 
divergence appears. The cone has two normal vectors in M}'^, given by = (1,0,0,0) and 

= (0, 0, r, 0) in {t, r, 6, (p) coordinates. Only the second fundamental form associated with 



is non- vanishing, with non-zero component fc^. = ^rsm2Q. The c-anomaly term in (1.1) 



then contributes 



c 



4807r 



3 (5.21) 



C COS^ fl , 2 



240 sin(] 



log^ € + ■■ 



In going from the first to the second line in (5.21), we assume that the tip of the cone is 



cut-off at some short distance tq oc e. Note that the a-anomaly does not give an additional 



contribution to the singularity because J = 0. The new UV divergent term (5.21 ) vanishes 
at f2 = |, where there is no conical singularity, while its coefficient diverges as Q goes to 
zero. 

A more heuristic way to obtain the log^ e term, similar to an argument for the cusp 
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geometry in p] is as follows. When f2 is small, the cone may be approximately decomposed 
into a union of cylinders with radius R = LsinQ and the length AL, where L is the length 
from the apex of the cone to one of the cylinders and AL is supposed to be small. The 
logarithmic term of the entanglement entropy of the cylinder comes from the c-anomaly 



given in (1.1) 



r AT 

5cynnder=^-^l0g(i?/6). (5.22) 

It follows that the entanglement entropy of the cone has the square of the logarithmic 
divergence 



r^max C 1 C 1 



which reproduces the leading behavior of (5.21) in the small O limit, except it is smaller by 



a factor of 2. We will see below that the factor in (5.23) is correct. 

We now present a more precise holographic derivation of the log^ e divergence, which 
correctly takes into account the conical singularity. We use the AdS^ metric 

o dy'^ - dt^ + dr'^ + r'^idO'^ + sin^ Odd)) ,^ 
ds^ = — — , (5.24) 

y 

so that the entangling surface is given by the cone S = {0 < r < oo, 6 = Q, < (p < 27t} at 
the AdS boundary y = 0. Taking the holographic coordinate y to be a function of (r, 6*), we 
use the conformal and rotational symmetries of AdS spacetime and the entangling surface 
to restrict the ansatz to 

l/(r,^) = ^. (5.25) 
The entanglement entropy functional is then given by 

^=-^) j- I d6 sin 9 ~gy^^^TF+W 

4:G AT J ^ J 



dgg^{g^ + g^)s'{gY + l-s{gf 



(5) 

1 

IT 



2G^ 



2G^ 



(5.26) 



.(5) 
'N 

where we introduced the new variable s = cos 9, which runs from sq = cosf2 to unity, and 
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redefined g{s) = g{6{s)). In the last equality, we changed the integration variable from s to 
9- 

The entanglement entropy is given by evaluating the EE functional in the last line 



of (5.26) on the function s{g) which solves the Euler-Lagrange equation, 



Vi9' + 9')s'{gy + l-s{gy 



d 
dg 



9{9' + 9')s'{g) 



Vi9' + 9')s'{9y + l-s{gy 



(5.27) 



subject to the boundary condition s{g = r/e) = Sq, where the UV cutoff is put at ?/ = e in 
the AdS spacetime. However, to find the leading divergences in the entanglement entropy, 
it only is necessary to know the function s{g) near the boundary at g = r/e, which we 



may assume is a large number. Taking the large g limit of (5.27), one may verify that the 
asymptotic expansion for s{g) near the boundary is 



s{9) = ^o + j\ + 



9' 



(5.28) 



While the solution above only satisfies the boundary conditions up to a term of order (e/r)"^, 
the difference does not affect the leading two singular terms in the EE. In evaluating the 



entanglement entropy functional in (5.26) on the solution s{g) in (5.28), we must evaluate 
the integral 



50 



dg 



r/e 



^-49 



l~sl9 



l^^Alog^ 



+ ^ ■ ^ log(r/e) + 0(e°) 



8 sinf2 



(5.29) 



where go is the minimum value of g{s). Then, preforming the r integral from r = g^t to 



we obtain the entanglement entropy for the cone: 



5, 



cone ^t;\ 



Ay, IT COS^ Q 



2e2 8sin(] 



log^(rmax/e) + 



(5.30) 



In order to compare this result with the naive calculation in (5.21), we use (5.13). One 



then sees that (5.30) is smaller than (5.21 ) by a factor of 2. As stressed above, the approach 



of (1.1) is not precise for singular entangling surfaces. It is nice, therefore, that it is only a 



factor of 2 off from the precise holographic result (5.30). 
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A Perturbed CFT 

In this appendix we discuss the entanglement entropy of (2 + l)-dimensional CFTs which 
have gravitational duals, perturbed by relevant operators. Many of these results were recently 
found in (24|. 

We work with an effective (3 + 1) -dimensional gravitational theory, with metric 

dsl = % f ^ - dt' + dr' + r'deA , y>0, (A.l) 

r \f{y) J 

and we assume that in the UV (small y) the metric asymptotes to AdS^ with radius L^y, 
i.e. f{y) = 1 -|- 0{y°')^ a > 0. The entanglement entropy across a circle of radius R in the 
boundary QFT is then given by the area functional 

where the function r{y) satisfies the Euler-Lagrange equation, e is the short-distance cut-off, 
and is the maximal value of y for the solution. 

If the effective (3-|-l)-dimensional gravitational theory comes from an exact D-dimensional 



string or M-theory background, with metric as in (3.1), then we may identify 



Let us begin with the conformal limit, where we may take f{y) = 1. We also define the 
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(3 + l)-dimensional Newton's constant G^^ = CjvfV^uv, where, in 10- and 11-dimensions, 
the Newton's constant takes the values G^^^^ = Svr^a'^^f^ and G^^^^ = (327r^)^^(27r£p)^, re- 
spectively. The solution to the equation of motion is then r{y) = W' — y"^. Evaluating the 
area functional on this solution and expanding in e gives 



S{R) 



ttL 



uv 



2G 



R 



R 



y 



2G^^e 



2G 



(4) • 

N 



(A.4) 



Suppose that the 3-dimensional CFT comes from the near horizon limit of M2-branes 
at the tip of the cone Cy = M x F, where Y is some 7-dimensional internal Sasaki-Einstein 
space. In the large A^ limit the theory is well described by the supergravity background 



2 r/s^ 



^UV 



V0Ud54 



*nF4 = 384L^vVoly , (A.5) 



where the radius Luv of AdS^ is quantized in plank units: 



A^ 



(2vr£p) 



6Vol(F) L 



uv 



£6 



(A.6) 



Substituting the relation between A^ and Lyy in (A.6) into (A.4) gives the renormalized 
entanglement entropy 



uv 



27r6 



^' V27Vol(r) 



+ 0(Ari/2) 



(A.7) 



As a consistency check, we may verify that this is equal to the finite part of the Euclidean 



free energy of the theory on the 3-sphere 21 ,47 , in agreement with the general result of 19 



Now suppose that there is an RG flow in the boundary QFT caused by perturbing the 
UV CFT by a relevant scalar operator O of dimension 1/2 < A < 3. The operator O is dual 
to a massive scalar field in the bulk, with A(A — 3) = m^. The bulk action is then given by 



1 



167rG 



(4) 
N 



d X 



-9 



R + 



6 



1 



-^uv 



1 



2 i2 

m (p 



where the dots stand for higher order terms in and the contributions of other fields, which 
won't be relevant for this discussion. When the field vanishes the equation of motion for 
the metric simply gives AdS^. The equations of motion for in the AdS^ background gives 
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the following asymptotic solution at small y, 



ct){y,x)=y^-^[Mx) + 0{y')], 



(A.9) 



where we have chosen the solution which stays finite at the boundary. This is the solution 
which corresponds to perturbing the action of the UV boundary CFT to 



(A.IO) 



The field (j) has a back-reaction on the metric. In particular, when (j)o{x) 
we find that at small y the Einstein equation gives us 



f{y) 



+ ■■■ 



is constant, 



(A.ii: 



To find the first correction to the renormalized entanglement entropy as a result of the 



relevant deformation of the UV CFT, it is sufficient to evaluate the action in (A. 2) with f(y) 
given in (A. 11) on the UV solution r(?/) = \J B? — y"^. A straightforward calculation then 
gives 



(A. 12) 



The factor 3 — A plays an important role; it ensures that the renormalized entanglement 
entropy is not changed by marginal perturbations. 
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